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Abstract - A difference analogue of the logistic equation, which preserves integrability, is 
derived from Hirota's bilinear difference equation. The integrability of the map is shown 
to result from the large symmetry associated with the Backlund transformation of the KP 
hierarchy. We introduce a scheme which interpolates between this map and the standard 
logistic map and enables us to study integrable and nonintegrable systems on an equal 
basis. In particular we study the behavior of Julia set at the point where the nonintegrable 
map passes to the integrable map. 



1. INTRODUCTION 

The completely integrable systems have been revealed to play a fun- 
damental role in the study of mathematical physics[l,2]. They provide 
a common basis to understand phenomena in various fields of mathe- 
matical physics. If we try to deform the systems but preserving their 
integrability we must generalize the symmetry possessed by the systems. 
The quantum group is such an example [3]. The essential feature of the 
new symmetry is to replace the differential operators by their difference 
analogue [4]. On the other hand it has been known [5] that many of the 
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completely integrable systems have their discrete counterparts which are 
also integrable. A typical example of this sort is the Toda lattice which 
is a discrete analogue of the KdV equation along the spatial direction[6]. 

It is, however, not obvious that equations obtained by making dis- 
cretization of variables of an integrable differential equation are inte- 
grable. The KdV equation can be derived from the Toda lattice by taking 
an appropriate continuous limit of the spatial variable. But, if we start 
from the KdV equation and discretize the spatial variable we obtain a 
variety of difference equations besides the Toda lattice. Those equations 
obtained by discretizing the variables in different ways do not preserve 
their integrability but create chaos in general. Our main concern in this 
article is to characterize the discretization of variables in integrable sys- 
tems, which maintain their integrability. 

Starting from the (2+1) dimensional Toda lattice and discretizing the 
two time variables such that all three variables in the equation appear 
symmetric, we arrive at Hirota's bilinear difference equation (HBDE)[7]. 
This equation is known completely integrable and is equivalent to the 
infinite number of soliton equations belonging to the KP hierarchy [8]. 
In other words, Hirota's difference equation is a simultaneous discrete 
counterpart of soliton equations in the KP hierarchy. As a matter of fact 
we can derive all of the soliton equations in the KP hierarchy from HBDE 
by taking appropriate continuous limits of the variables. Therefore HBDE 
is a desirable object to study for our purpose. 

A system which evolves along discrete time can be considered as a 
dynamical system if we regard it as an infinite sequence of maps. In this 
sense our problem might be thought equivalent to characterize dynamical 
systems which are integrable. Since the dependent variable of HBDE is 
defined on an arbitrary Riemann surface, we are especially interested in 
the complex dynamical systems. 

From physical point of view the existence of discretization of inte- 
grable systems which preserve integrability is important by itself. The 
intimate correspondence between the solvable statistical models on two 
dimensional lattice and the two dimensional conformal field theories has 
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been noticed from the begining of the study of the integrable systems [9]. 
In the relativistic quantum field theory there arises infinity owing to the 
fact that the space-time is a continuous media of the fields. The lattice 
gauge theory[10] provides a way to avoid this difficulty The space-time 
is replaced by a set of lattice points. Only at the end of all calculation of 
physical quantities the lattice constant is set zero. This approach causes 
many problems specific to the way the lattice constant is introduced into 
the theory. The relativistic invariance, for instance, is lost in general. If 
the theory on the lattice is integrable, however, we are certain that the 
results obtained on the lattice are correct in the continuum limit. The 
systems which preserve integrability after discretization of variables are 
also peculiar in the sense that all results calculated on computer are exact 
without errors. 

In order to clarify the notion of integrability preserving discretization 
of integrable systems we like to study, in this paper, the simplest case of 
Hirota's bilinear difference equation. We first show that if we consider 
the dependence on one variable and ignore the other variables, HBDE 
reduces to a difference analogue of the logistic equation, which preserves 
integrability. Secondly we argue that the discrete time evolution of this 
integrable map can be considered as a special case of the Backlund trans- 
formation. The existence of Backlund transformation characterizes inte- 
grable systems [11]. It enables one to derive a sequence of new solutions 
starting from a known solution. Hence it is quite natural to intepret that 
the Backlund transformation generates a time evolution of the integrable 
system. In the last part of this paper we introduce a scheme which en- 
ables us to discuss the completely integrable systems and nonintegrable 
systems on the same ground. Specifically we consider analytic properties 
of a map which interpolates two different types of discretization of the 
logistic equation. Namely we investigate a map which combines one of 
the standard logistic map, which presents chaos, and the one derived from 
HBDE through the reduction, which is completely integrable. 
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2. DISCRETIZATION OF LOGISTIC EQUATION 



We are interested in characterizing integrable systems which preserve 
their integrability after discretization of dependent variables. Most of 
the integrable systems discussed in the literature of theoretical physics 
have infinite number of freedom and thus too much complicated to be 
considered as dynamical systems. Instead of studying them directly we 
will start our discussion in this paper from the simplest case, i.e., one 
dimensional ordinary differential equations. One of such equations is 
given by 

d i - '<*>■ m 

This is equivalent to one dimensional Newton equation, which can be 
solved as 



There exists a unique solution as long as l/f(x) is integrable. If we dis- 
cretize the variable t according to Euler's prescription, we obtain, writing 
x(t) -> Zf. 



When h — 1, this equation turns to be an infinite sequence of a map, 
i.e., a one dimensional dynamical system, h — > reproduces the original 
equation. If we consider z% on the complex plane, and f(z) is a rational 
function of z with degree larger than 2, this map is known[12] to have 
a Julia set (the closure of repelling periodic points). Therefore this map 
shows chaotic behaviour somewhere on the complex plane, and is not 
integrable. When / is a polynomial of z of degree 2, the map is called a 
logistic map and is studied in detail[12]. In fact if we put 




(2) 



zi+h = zi + hf(zi). 



(3) 



f(z) = z(ii-h- jiz) 



(4) 



and h — 1 we get the standard form of the logistic map: 



Z l+1 = fJ,Zi(l - z t ). 



(5) 
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In the limit of h — > it turns to the logistic equation: 

^=fix(l-x). (6) 

The discretization of the logistic equation is not unique. Hence there 
arises a question; what is the integrable discretization of the logistic equa- 
tion? To answer this question let us recall some properties possessed by 
integrable discrete systems. Hirota's bilinear difference equation (HBDE) 
is one of such equations whose solutions are completely known[8,ll]. It 
is given by 

afiP + 1, Q, r)f(p, q + 1, r + 1) + (3f(p, q + 1, r)f(p + 1, q, r + 1) 

+ I f(p,q,r + l)f(p+l,q + l,r) = 0. (7) 

Here p, q, r are discrete variables and a, ft, 7 are arbitrary parameters 
subject to the constraint a + (3 + 7 = 0. Infinitely many equations of the 
KP hierarchy can be reduced from this single equation by taking continu- 
ous limits of various combinations of the variables and the parameters [7]. 
To be more specific we introduce a new set of variables {t n } = {ti,t 2 ,- ■ ■} 
by[8] 

tn = -T, k J Z h U ^ 1 - ( 8 ) 

n j 

Here we assume Y,j kj = 0, which is nothing but the energy- momentum 
conservation law in the language of the string model[2]. This new set of 
variables {t n } are the independent variables of the soliton equations in 
the KP hierarchy. The variables p, q, r of HBDE are any three of k/s of 
(8). Therefore we consider infinite number of copies of HBDE. To obtain 
a particular soliton equation we change, in one of the HBDE, the variables 
from {kj} to {t n } and reduce the number of t variables by neglecting their 
dependence. From this argument it is apparent that a large symmetry of 
the KP hierarchy can be seen if we use the variables k/s instead of t n 's 
of the soliton equations. 

The solutions of HBDE are known explicitly and given by analytic 
functions defined on Riemann surface with arbitrary genus[13]. {zj} are 
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local coordinates of the Riemann surface. The space of the solutions is 
called universal Grassmannian which is a Grassmann manifold of infinite 
dimension [14]. The equation which characterizes such a large space of 
solution is rather simple if we consider HBDE instead of the set of soliton 
equations in the KP hierarchy. 

If we wish to discuss integrable systems in comparison with noninte- 
grable systems such as logistic map, however, HBDE is still complicated. 
Therefore we try to reduce the dependence on the variables. As we sim- 
ply neglect dependence on two independent variables, say q and r, HBDE 
becomes trivial. In order to obtain the most simple nontrivial reduction 
we introduce new set of variables 

3 1 
l = p + q + r + -, m = -q--, n = p + q+l (9) 

and denote by g the dependent variable written by new variables as 

f(p,q,r) = g n (l,m). 

We then obtain 

ag n (l + l,m)g n (l, m + 1) + f3g n (l, m)g n (l + 1, m + 1) 

+ 70 n+ i(Z + l,m)0 n _i(Z,m + l) = O. (10) 
The first nontrivial relation is given if we choose 

9l (l,l) = g (l,2) = <? 2 (U) = zi 

0i(Z,2) = c(const.) (11) 
g n {l,m) = 0, otherwise 
and also (3 = —fxa, 7 = (/i — l)a, c = ^— ^, from which follows 

zi+i = /izi(l - zi+i). (12) 

This is quite similar to the standard logistic map (5) except for the second 
term on the right hand side. Since this equation is derived by a reduction 
starting from an integrable difference equation, the behaviour of solutions 
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must be completely different from those of the standard map. In fact 
writing (12) as 

z l+l = T— 13 

1 + fXZi 

the solution is given by 

Zl = l + u^z (14) 

for arbitrary initial value z . This solution is exactly the one which co- 
incides with one of the logistic equation (6) in the continuous limit. It is 
easy to check that this equation itself reduces to the logistic equation in 
the same limit. In this sense we call (12) the integrable discretization of 
the logistic equation. 

The equation (12) itself was studied by Morisita[15] as a difference 
analogue of the logistic map and revived by Yamaguti[16] and Hirota[17]. 
We like to emphasize here that the integrability of the map (12) has the 
same origin as the integrability of the KP hierarchy. 



3. BACKLUND TRANSFORMATION 

As we have seen in the previous section a slight difference of the way 
of discretization of an integrable equation leads to a big difference of the 
behaviour of the system. The question we ask in this section is how a 
regular map is possible and does not fall into chaos. To this end we will 
show that the map which maintains regularity is a particular sequence 
of the Backlund transformations of the KP hierarchy. In the case of the 
KP hierarchy the Backlund transformations are auto Backlund transfor- 
mations[7,18]. In other words the transformation generates a sequence 
of solutions of the same equation. Starting from a point of the universal 
Grassmannian a sequence of the Backlund transformations sweeps a tra- 
jectory on the space of solutions. From this result we can glance the deep 
connection of the integrable dynamical systems and the large symmetry 
associated with the Backlund transformation. 

To begin with we first recall the fact that the time variable / of the 
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map (13) is a certain combination of k/s. We like to show that the 
Backlund transformations of the KP hierarchy are the operations which 
change values of k/s by one. It then follows that there exists a certain 
transformation which shifts / into I + 1 without changing other variables. 

The Backlund transformation of the KP hierarchy is generated by the 
operator [11] 

B(z,w)=A(z)A(w) (15) 



where 



OO OO 1 

A(*) = exp[- £ t n z~ n ] exp[£ -z n d n ], (16) 



n=l n=l 



OO OO 1 

A(z) = expE t n z~ n } exp[- £ -z n d n \. (17) 

n=l n=l U 

In this expression d n means d/dt n , where t n 's are those variables of the 
soliton equations in the KP hierarchy. Now we consider the operation of 
A(zj) to f(ki, /c2, • • •)• The right hand side factor of A(zj) transforms t n 
according to 

tn -> *n+ -2?, Tl > 1. 

n J 

We observe from the expression of t„ (8) that this change of variables 
amounts to replace kj by kj + 1. Similarly the left hand side factor of 
A(zj) can be rewritten in terms of k/s. In doing this we must be a bit 
careful because A(z) is not well defined when z approaches to one of z/s. 
Thus we define a new regularized operator by 

\{*i) = U(*i (18) 
i 

A(zj) is recovered in the q — > 1 limit. From this expression we see that 
the operation of B q (z i: Zj) = A q (zi)A q (zj) to one of solutions of HBDE, 
say /(fci, fc 2 , • • •), yields 

B q (zi, Zj)f(k\, ■ • ■ , ki, ■ • • , kj, • ■ ■) = 6(z«, Zj)f(ki, • • • , • • • , kj—1, • • •), 

(19) 

where 

■Zi-qzi\h 1 



^•)=np^) 



J — 1 Z l Z j ~ Q Z i 
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Therefore the operation of B of the Backlund transformation (15) causes 
to f(ki, k 2 , ■ ■ •) the shifts of the variables ki — * k{ + 1, kj — > kj — 1. 
Let f^(ki, k 2 , • • •) be the function obtained from /(fci, &2, '' ') by the 
Backlund transformation, i.e., 

e aB q ( Zi , Zj)f{k ^ ^ . . , } = /( a) (A;i> fca> . . (2Q) 

It is well known that the generator of the Backlund transformation i? has 
a fermionic nature, i.e., B 2 (zi, Zj) = holds [11]. Hence we can write 

aB q ( Zi , Zj)f{k u k 2 , ■■■) = f^ih, k 2 , ■■■)- f(k ± , k 2 , ■■■). 

To be more specific let us choose four out of {kj} and call them 
p,q,r,s. Further we define l,m,n from them according to (9). Then 
we see that the increase of I, m, n by one are equivalent to the change of 
p, q, r as 

I -> / + 1 <^> r -> r + 1 

m— >m+l44>p— q q — 1 (21) 

n — > n + 1 4=> p — > p + 1, r — > r — 1. 

If s is the variable associated to the singular point z s on the Riemann 
surface, the amplitude / does not depend on s. Letting z/s associated 
to p, q, r, s as z p , z q , z r , z s , the shifts of (21) are realized by the following 
operations: 

I -> I + 1 44> B 3 (z r , z 8 ) 

m m + I Bg(z p , z q ) (22) 

7i — > n + 1 B q (z p , z r ). 

From this we see that under the operation of B q (z r , z s ) only the variable 
/ changes its value whereas other variables remain unchanged. According 
to (20) we can express g n (l + l,m) by gj^(l,m) and g n (l,m). In other 
words g n at Z+l is determined by and g n ed,l. Hence it is the generator 
of the Backlund transformation which generates the time sequence of zi 
under the integrable map (13). 
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Finally we notice about the symmetry corresponding to the Backlund 
transformation. If we expand B(z,w) into double Laurent series as 

B(z,w)= J2 B mn z m w~ n (23) 

m,n£Z 

we can show [11] 

\B mri) B m i n i\ 5 m ' n B mn i 5 mn iB m i n . (24) 
Hence B generates a group GL(oo). 



4. AN INTERPOLATION OF INTEGRABLE AND NONINTEGRABLE 
DIFFERENCE ANALOGUES OF THE LOGISTIC EQUATION 

In the previous sections we derived the integrable map associated with 
logistic equation starting from HBDE and showed that the integrability 
of the map owes to the symmetry associated to the Backlund transforma- 
tion of the KP hierarchy. We like to clarify in this section how integrable 
maps can be characterized among nonintegrable ones. For this purpose 
we study a map which interpolates between the integrable and the non- 
integrable difference analogues of the logistic equation, and argue the 
behaviour of the Julia set. 

We have seen in §2 that we can derive two different maps (5) and (12) 
starting from the same logistic equation (6). Although they look similar 
at the level of eqation the behaviour of their solutions are quite different. 
We like to know analytically the transition between these two systems. 
To proceed we generalize (5) and (12) such that it is a discretization of 
the logistic equation (6) and (5) and (12) are included as two particular 
cases. A simple generalization is given by 

zj+i =fizi{l -'yzi - (1 -7)z,+i}, (25) 

where 7 is an arbitrary parameter. It is apparent that the standard 
logistic map and the integrable map are recovered if 7 is set 1 and 
respectively. If we solve this equation for zi+i we find 

zi+i = <l>(zi), (26) 
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where we have introduced a map (f>(z); 



(27) 



l + fl (l-j)z 

This is what we are going to analyze in this section. 

We first notice that the map 4>(z) of (27) is similar to the one known 
as Blaschke product of degree two [19]; 

e *3±t± (28) 
1 + Az' K J 

where 8 is an arbitrary real number and A is the complex conjugate of A. 
This maps the unit circle on itself in the complex plane. Moreover the 
Julia set of this map is the unit circle. Since there is no other unstable 
fixed point in this map it is very convenient to concentrait our attention 
to this to analyze (27) for the first step. 

Comparing our map with the expression of the Blaschke product, we 
see that the map (27) turns to (28), a Blaschke product, if the parameters 
satisfy the relations; 

„ 7 =-e», ^-^y (29) 

We are interested in the case of small values of |7|. As far as (29) are 
satisfied, the Julia set remains on the unit circle. This happens only if 
the value of fj, diverges when 7 approaches to zero. Hence we can not 
draw any conclusion for the limit of 7 — > 0, since the value of \i in our 
interest is finite. 

In order to proceed further we try to see the map from different coor- 
dinate system in the complex plane. Namely we introduce a new variable 
w which is a Mobius transformation of z, i.e., w = ip(z). In this new 
coordinate system, <f> turns to following to 

<f)(z) -> <f)(z) =ip- 1 o(f)o^(z). (30) 

The map (f>(z) of (27) has three fixed points at 0, 00 and 1 — K We want 
and 00 to remain as fixed points after the change of the variable, so 
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that the new map still keeps the form 



^) =e ^J^+£_. (31) 



One of such coordinate transformations is 



1-i 

M z ) = ^_ , (32) 

which yields a new map of the form (31) with 

A= l + 7^7 e _* y = A , e * (33) 
/i + 7 - /i 7 

In the new coordinate the fixed points are at 0, 00 and 

1 



V = 



+ 7 - /ry 

The multipliers of the fixed points at 0, 00 and p are given, respectively, 

as 

\ 1+7-/^7 , XI 1 fOA\ 

A = — ■ , Aoo = /i, Ap = 1 . (34) 

fi + 7 - m 7 

Notice that z = p is a repelling fixed point irrespective to the value of 

/i when I7I is small. The nature of the fixed point at 00 is uniquely 

determined by the value of \x alone. 

If /i varies independent on the value of 7, the product Ao • Aoo goes to 

1 as 7 tends to 0. This means that, for small values of I7I, the origin is 

attractive when 00 is repelling and vice versa. This situation can be also 

understood from the map itself. Letting 7 be zero in (31) we get 

ij>(z) = -z. (35) 

The sequence of this map attracts all points in the complex plane to the 
origin (resp. 00) if is larger (resp. smaller) than 1. Therefore the 
Julia set reduces to a point at 00 or depending on the value of is 
larger or smaller than 1, respectively. 

The condition for the map (31) to be a Blaschke product is 

1 + 7-/^7 - M 2 - 1 

— . =n or 7 = — . 36 

/i + 7-^7 \ii - l 2 
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Under the constraint (36) for the parameters we calculate the values of 
the multipliers (34). We then find Ao = p, \<x> = P, hence the origin as 
well as the oo are attractive fixed points when \p\ < 1. The third fixed 
point is at 



Therefore p is on the unit circle and is a repelling point where 7 takes 
a negative real value. It then follows that all points inside of the unit 
circle are attracted to the origin by the map whereas the points outside 
of the circle are pushed away to infinity. The behaviour of the map on 
the unit circle is chaotic, because the Julia set is the unit circle. It can 
not escape from the unit circle as long as p satisfies the constraint (36) 
for a given value of 7. This is true also when the value of I7I approaches 
to zero. Does this mean that the Julia set remains at 7 = ? To see this 
point more carefully, we restrict the value of p to satisfy (36) and write 
the map explicitly; 



Now let 7 go to zero along the negative real axis, so that \p\ < 1 is 
satisfied. The absolute value of p tends to one owing to the constraint 
(36). The multipliers at and 00 approach to 1, and thus these points 
cease to be attractors. Then the map is given by (35) again as we see 
from (37). In this case, however, p is pure phase and the sequence of this 
map forces every point on the complex plane to rotate uniformly around 
the origin. The angle of the rotation is determined by the phase of p. 
From (36) we can calculate the value of the phase. Writing p = \p\e t27TX 
it turns out to be 



Therefore the manner of the rotation is uniquely determined by the value 
of the limit of the ratio of 1 — \p\ and 7. 

It is known that a sequence of the rotation of an irrational number 
of phase causes one dimensional ergodic behaviour of zero entropy [20]. If 




(37) 



2ttx = lim arccosf 1 + 
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we choose the ratio between 1 — and 7 arbitrarily the phase of \i will 
fall to an irrational number with the probability one as 7 approaches to 
zero. Then the rotation of all z on the complex plane caused by (37) are 
most probably ergodic with entropy 0, except for the neutral fixed points 
at and 00. The Julia set, which was on the unit circle when 7 < 0, does 
not exist when I7I is set equal zero. 

Let us summarize the above results. As the value of 7 approaches to 
we have three cases: 

1) if > 1, the Julia set is at 00, 

2) if < 1, the Julia set is at 0, 

3) if \fi\ = 1, there exists no Julia set but all the points between and 

00 on the complex plane rotate ergodically so long as the ratio between 

1 — \fj,\ and 7 does not go to a rational number as 7 goes to 0. 

We can further proceed the analysis of our map (27) from the math- 
ematical point of view. As a complex dynamical system the map is in- 
teresting by itself. Since further investigation requires more preparation, 
however, we shall report it elsewhere [21]. 

DISCUSSIONS 

The completely integrable systems constitute a small portion excep- 
tional among general nonlinear systems and behave quite different from 
others. Nevertheless they describe important nature of fundamental prob- 
lems in various fields of theoretical physics. The large symmetry of the 
integrable systems enables us to determine their solutions analytically, 
whereas the missing of some conserved currents prevents us to predict the 
behaviour of solutions in the generic case. The soliton type of equations 
are completely integrable systems whose detail features have been studied 
in the last two decades. In order to proceed further and to understand 
the real importance of the role played by integrable systems in theoretical 
physics it is necessary to study them in comparison with nonintegrable 
systems. In other words the real understanding of the basic problems in 
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theoretical physics will be achieved when the complete integrable systems 
are fully characterized among general nonlinear systems. 

There has been, however, no systematic way to investigate noninte- 
grable systems. The nonintegrable systems which have been studied are 
mainly one dimensional dynamical systems. Even the simplest one of such 
systems is far from totally investigated. This should be contrasted with 
the case of the completely integrable systems, in which the most interest- 
ing systems are infinite dimensional and solutions are given analytically. 
Therefore it is difficult to find a way to discuss these two different systems 
on a common basis. 

In this paper we have shown that the simplest reduction of Hirota's 
bilinear difference equation turns out to be an integrable difference ana- 
logue of the logistic equation. We provided a scheme which interpolates 
this equation and the standard logistic map. It enables us to study in- 
tegrable and nonintegrable maps on the equal ground. In particular we 
concentrated our attention to the behaviour of the Julia set at the point 
where the system turns from nonintegrable to integrable. 

The integrable part we studied here, is a small portion of HBDE. It 
will be the most desirable if there is a simple scheme in which HBDE itself 
is included as an integrable part of some nonlinear system. Although we 
have not found such a scheme we could interpret the integrability of the 
map results from the large symmetry possessed by the integrable system 
such as GL(oo) which generates the Backlund transformation of the KP 
hierarchy. 



15 



REFERENCES 



1. See, for example, "Infinite Analysis", ed. A.Tsuchiya, T.Eguchi and M.Jimbo 
(World Scientific, Singapore, 1992); "Yang-Baxter Equation in Integrable Sys- 
tems", ed. M.Jimbo (World Scientific, Singapore, 1990). 

2. S.Saito, Phys. Rev. D36 (1987) 1819; Phys. Rev. Lett. 59 (1987) 1798; Phys. 
Rev. D37 (1988) 990; H.Cateau and S.Saito, Phys. Rev. Lett. 65 (1990) 
2487; in "Strings '88" ed. S.J.Gates,Jr., C.R.Preitschopf and W.Siegel (World 
Scientific, 1989) p. 436 ; in "Nonlinear Fields; Classical, Random, Semiclassical" , 
ed. P.Garbaczewski and Z.Popowicz, (World Scientific, 1991) p. 268. 

3. M.Jimbo, Lett. Math. Phys. 10, (1985) 63; 11, (1986) 247; Int. Jour. Mod. 
Phys. A4 (1989) 3759; V.G.Drin'fcld, Sov. Math. Dokl. 32 (1985) 254; 
"Quantum Groups", ICM Proceedings, (Berkeley, 1986) P.898. 

4. I.B.Frcnkcl and N.Yu.Reshetikhin, Comm. Math. Phys. 146 (1992) 1. 

5. E.Datc, M.Jimbo and T.Miwa, J. Phys. Soc. Jpn. 51 (1982) 4116, 4125. 

6. M.Toda, J. Phys. Soc. Jpn. 22 (1967) 431; "Theory of Nonlinear Lattices" 
(Springer-Verlag, 1981); N.Saitoh, J. Phys. Soc. Jpn. 49 (1980) 409. 

7. R.Hirota, J. Phys. Soc. Jpn. 50 (1981) 3787; "Nonlinear Integrable Systems" 
ed. M.Jimbo and T.Miwa (World Scientific, 1983) p.17. 

8. T. Miwa, Proc. Jpn. Acad. 58A (1982) 9. 

9. A.Belavin, A.Polyakov and A.Zamolodchikov, Nucl. Phys. B241 (1984) 333. 

10. See, for example, C.Rebbi, "Lattice Gauge Theories abd Monte Carlo Simula- 
sions", (World Scientific, Singapore 1893). 

11. E.Date, M.Jimbo, M.Kashiwara and T.Miwa, J. Phys. Soc. Jpn, 50 (1980) 
3806, 3813; M.Jimbo and T.Miwa, Publ. RIMS. Kyoto Univ. 19 (1983) 943. 

12. See, for example, R.L.Devaney, "An Introduction to Chaotic Dynamical 
Systems", 2nd ed. (Addison- Wesley, 1989). 

13. I.M.Krichever, Russian Math. Surveys 32 (1977) 185; 
M.Mulase, J. Diff. Gcom. 19 (1984) 403. 

14. M.Sato and Y.Sato, Publ. RIMS. (Kyoto Univ.) 388 (1980) 183 ; 
ibid. 414 (1981) 181; M.Sato, ibid. 433 (1981) 30. 

15. M.Morisita, Res. Popul. Ecol. VII (1965) 52. 

16. MYamaguti, Mathematics on Nonlinear Phenomena (in Japanese) 
(Asakura-Shoten, Tokyo 1972). 



16 



17. R.Hirota, Applied Mathematics 3 (in Japanese) (1993) 48. 



18. S.Saito and N.Saitoh, J. Math. Phys. 28 (1987) 1052; Phys. Let. A 123 
(1987) 283; N.Saitoh and S.Saito, J. Phys. Soc. Jpn 56 (1987) 1664; J. Phys. 
A:Math. Gen. 23 (1990) 3017. 

19. See, for example, J.Milnor, "Dynamics in One Complex Variable: Introductory 
Lectures" SUNY preprint #1990/5 (1990). 

20. P.Billingsley, "Ergodic Theory and Information" (John Wiley and Sons. Inc., 
New York, 1965). 

21. N.Saitoh, S.Saito and A.Shimizu, "An Analysis of a Family of Rational Maps 
which Contains Intcgrable and Nonintegrable Difference Analogue of the Lo- 
gistic Equation" , in preparation. 



17 



